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Abstract 

We propose a synthetical weights' dynamic mechanism for weighted networks which 
takes into account the influences of strengths of nodes, weights of hnks and incom- 
ing new vertices. Strength/ Weight preferential strategies are used in these weights' 
dynamic mechanisms, which depict the evolving strategies of many real- world net- 
works. We give insight analysis to the synthetical weights' dynamic mechanism and 
study how individual weights' dynamic strategies interact and cooperate with each 
other in the networks' evolving process. Power-law distributions of strength, degree 
and weight, nontrivial strength-degree correlation, clustering coefficients and assor- 
tativeness are found in the model with tunable parameters representing each model. 
Several homogenous functionalities of these independent weights' dynamic strategy 
are generalized and their synergy are studied. 
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1 Introduction 



Complex networks pipilHllHIS] depict a great many real- world networks like the 
scientific collaboration networks (SCN) P71l8]l9] . World Wide Web (WWW) [10], 
world-wide airport networks (WAN) [TT][T2] and so on. Simple binary net- 
works pl2|3] are used to depict the topological aspects of these real-world 
networks. Degree distribution and degree related clustering coefficients can be 
analyzed from the model. Typically, Barabasi and Albert proposed a linear 
preferential attachment model on which most other binary models are based 
on (BA model[T3]). However, real- world networks often contains far more in- 
formation than that binary networks can express, since relations in the net- 
works are not necessarily binary. Therefore links with weights are introduced 
to emphasize the importance of heterogenous relations between nodes in net- 
works. Barrat, Barthelemy and Vespignani first build a model for weighted 
networks based on preferential attachment mechanism (BBV model [T^fTS] ). 
Since then various models have been derived to mimic diverse real-world net- 
works |^|f5|16|17|18|19|20ll21|22] . all of which emphasize a different kind of 
weights' dynamic mechanism. 

The flourishing research on various weights' dynamics mechanisms roots in the 
the heterogenous behaviors of real-life networks. Current research has already 
covered most part of typical behaviors of weights' dynamics. However, the 
interaction and operations of these scattered weights' dynamics are poorly 
studied and their underlying homogeneity are still not know. In our paper, we 
generalize weights' dynamics with monotonous weights' growth and focus on 
interactions and cooperations among them. 

With a close scrutiny into prevailing weights' dynamic models depicting real- 
world networks we find that there are three main sources of weights' increment 
dynamics: the variation of traffic caused by introducing of new nodes [HI5|16|17P1 
the increment of links' weights based on links' weights themselves [2^, and the 
increment of weights based on strength of two ends of the link |20f21ll22j . Most 
weights' increment dynamics can be grouped into these three sources. Many 
works have done to empirically validate the classification, among them New- 
man gives the most comprehensive and convincing experiments [23]. In this 
paper we will give numerical and analytical study to the synthetical weights' 
dynamic model which comprises the three above-mentioned mechanisms, and 
reveal the scale-free characteristics and nontrivial clustering coefficients and 
assortativeness of the model. 

Our paper is organized as follows. In Sec. 2 we define basic definitions and 
terms to represent the weighted network and gives a brief review of related 
works done. In Sec. 3 we depict our model and give clear definition of three 
weights' dynamic mechanisms. In Sec. 4 we analytically calculate the mathe- 
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matical expression for probability distributions of strength and weight, degree- 
strength correlations and related attributes. In Sec. 5 we preform simulations 
to mimic the proposed mechanisms and analyze the experiment results in 
detail. In Sec. 6 we give a conclusion to the paper. 



2 Definitions 

Weighted networks can be represented by a adjacent matrix W where Wij 
defines the weight of link between vertices i and j. Wij — indicates that 
there is no link between vertices i and j. Therefore topological and weighted 
information can both be revealed from W. Matrix W is symmetrical therefore 
Wij = Wji. Degree ki defines the number of vertices vertex i is linked with, and 
strength Sj defines the total weights of links that ends in the particular vertex 
i. Si can be written as Sj = J2j Wij, and ki can be written as ki = J2j sgn{wij), 
where sgn{) is the signum function. -P(s), P{k), P{w) defines the probability 
distribution of strength, degree and weight. Previous studies have revealed 
that in many weighted networks, P{s), as well as P{k) and P{w), displays 
a power-law distribution as P{s) = s~'^'' , P{k) = k~""' and P{w) = w^"*^ . 
There is also a power-law correlation between s and k that s = k^. Clustering 
coefficients and assortativeness of weighted network is also studied and the 
details will be discussed later. 



3 The Model 

The model proposed in this paper starts from an initial configuration of A^o 
vertices fully connected by links with weight wq = 1 (A^o-clique). At each time 
step, the network evolves under two coupled mechanisms: topological growth 
and weights' dynamics. Weights' dynamics are discussed in detail in this paper, 
where all three sources of weights' dynamics are taken into account. 

3.1 Topological Growth 

At each time step, a new vertex n is introduced into the network and con- 
nected to p existing vertices i. Vertices are chosen according to the strength 
preferential probability 

n„_, = , (1) 

and the weight of this new link is set to wq — 1- 
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3.2 Weights' Dynamics 



There are three sources of weights' dynamics: the local increment of weights 
triggered by the introduction of the new vertex, the self-increment of weights 
based on the weight of each link, and the mutual selection dynamics focusing 
on creation and reinforcement of links between existing vertices based on their 
strengths. These three weights' dynamics mechanisms interact and cooperate 
during the evolution of the network. There are several suggestive independent 
works for these three sources weights' dynamics: Barrat, Barthelemy, Vespig- 
nani first suggest the local rearrangement model considering the impact of in- 
coming vertices. Dorogovtsev, Mendes's work Wen-Xu Wang's works suggest 
the mutual selection model which well depict the second source. ********'g 
work initializes the idea of weights' self-increment although the idea is not 
comprehensively studied yet. There are lots of works done that empirically 
proving the validity of the three sources. Newman in his work by empirically 
studying the scientific collaboration network suggests that two scientists would 
have better chance to enforce their collaboration if they already have a lot of 
works together, and scientists are more likely to develop new collaborative re- 
lationships if they already have relatively large numbers of collaborators [24j . 
Collaborative relationships of scientists and their collaborators are analogous 
to weights and degree in a network. Degree can be generalized to strength 
if we take into account the amount of collaborations between each pairs of 
collaborators in stead of a binary expression. 

The introduction of new vertex brings variation in traffic across the network. 
For simplicity, we restrict the variation to the neighborhood of vertex i which 
has just been chosen to link with the new vertex. An overall increment of 
6 is introduction at each time step. The increment is distributed among the 
neighborhood of r{i) according to weight preferential mechanism: 

Wij Wij + 6— . (2) 

The strengths of i and all j G T{i) are also increased as a result of the increment 
of weights in the neighborhood of i. Considering the probability of vertices i 
been chosen, the increment of Wij can be rewritten as 

Awij = p— 5^ + p 6^ 

Z^k ^k Si 2^k ^k Sj 

= 2pS^ . (3) 



At each time step, n existing links are chosen to increase according to the 
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weight preferential probability: 

Wij Wij + n— (4) 

Each chosen link is increased hy wq — 1. The links with larger weights always 
have more chance to reinforcement. 

At each time step, each existing vertex i selects m vertices according to the 
strength preferential mechanism: 

There would be a alteration in links between i and j if and only if i and j 
have mutually selected each other. The probability that the linking condition 
between i and j changes can be defined to be: 

Sj 

n,; 1 = m— m 



J2k ^k ~ Si J2k ^k ~ Sj 

2 Si Sj 

m 



(EfcSfc - Si){J2k Sk - Sj) 

SiSj 



2 SiSj , , 

' . (6) 



If there is not a link between i and j, a new link with assigned weight Wq = 1 
will be added. If there is already a link between i and j, the link will be 
increased by Wq — 1. 

These three weights' dynamics mechanisms interact and cooperate during the 
the process of network development. Synthesize all the these three mecha- 
nisms, the increment of weights can be represented to be 



Noticing the fact that J2k,i ^ki = \ Efe Ski we can rewrite the above equation 
as 

Wij ^ Wij + {2p5 + 2n)-^ + m \'''' . (8) 

l^k Sk \l^k Sk) 



4 Evolution and Distribution of Degree, Strength and Weight 



Using the continuous approximation, we can assume that s, k, w, t are all 
continuous. Therefore we get 
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There are two sources contributing the increment of strength Sj, one is the 
weights' dynamic and the other is hnking with the new added node. Therefore 
the increment of Sj can be written as 



dsi X ^ dwij s 

-i7 = T.-jr + p 



dt ^ dt Efc-Sfc 



= {m^ + 2p6 + 2n + p)^^ . (10) 

The sum of strength of all nodes J^k^k at time t can be calculated as 

t 

Jo ^ dt 

= {m^ + 2pd + 2n + 2p)t , (11) 
and using this equation we can rewrite Eq. (ITUl) as 

dsi w? + 2p5 + 2n + p Si 



dt m? + 2p6 + 2n + 2p t 



(12) 



With the initial condition Si{t = i) = 1, we can integrate the above equation 
to obtain 

m'^+2pS + 2n+p 
ft \ m-^+2p5 + 2n + 2p 
Si{t) = \^-J (13) 

From the equation we can see that three parameters m, 5 and n cooperatively 
and interactively govern the growing speed of strength Sj. Is is really amazing 
to find out that all three sources of weights' dynamics influence the growing 
speed of strength in similar ways. The simulation of evolution of Si is given in 
Fig. [H We see how m, 6, n, p contribute to the evolution of Si independently by 
fixing three other parameters. We also show how these four parameters interact 
by varying them the same time as indicated by the above equation. We see Sj 
display a power-law distribution as t evolves, and variable m contribute larger 
alteration in Sj with relatively small amount of increment. 

The knowledge of the time evolution of the various quantities allows us to 
compute their statistical properties. The incoming time ti of each vertex i is 
uniformly distributed in [0,t] and the strength probability distribution can be 
written as 

1 



^(^'^) = 7-^ / ^(^ - ^iit))dU , (14) 
t + iVo Jo 

where d{x) is the Dirac delta function. Using Eq. f|T3|) we can get in the infinite 
size limit t ^ oo the distribution P{s) ~ s''^" with 

% = 2+ ^ . (15) 

m"^ + 2pd + 2n + p 
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Fig. 1. Evolution of degree with time for vertex i = 1, the time span is from 1 
to 5000 and the result is average of 10 individual experiments. Various values for p, 
m, n, 6 are chosen to display the interplay of different parameters. 
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Fig. 2. Evolution of degree ki with time for vertex i = 1, the time span is from 1 
to 5000 and the result is average of 10 individual experiments. Various values for p, 
m, n, 5 are chosen to display the interplay of different parameters. 

The degree probability distribution P{k) ~ k~'^'' can be obtained by combining 
s ^ k^ with Eq. (fT3l) . From the equation of the conservation of probabihty 



P{k)dk 



P{s)ds 



(16) 
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Fig. 3. Evolution of degree Wi with time for vertex i = 1 and j = 2, the time span 
is from 1 to 5000 and the result is average of 10 individual experiments. Various 
values for m, n, 6 are chosen to display the interplay of different parameters. 
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Fig. 4. Nontrivial correlation between degree k and strength s, the size of the network 
is 5000 and the result is average of 10 individual experiments. Various values for p, 
m, n, 6 are chosen to display the interplay of different parameters. 

we can get 

P{k) - P{s)^ ~ s-^^l3k^-^ ~ , . (17) 
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Fig. 5. Power-law probability distribution of strength s, the size of the network is 
5000 and the result is average of 10 individual experiments. Various values for p, 
m, n, 6 are chosen to display the interplay of different parameters. 




Fig. 6. Power-law probability distribution of degree k, the size of the network is 
5000 and the result is average of 10 individual experiments. Various values for p, 
m, n, 6 are chosen to display the interplay of different parameters. 

Therefore we get 7^ = /5(7s — 1) + 1 in P{k) ~ k~'^''. The simulations of ki and 
P(A;)are given in Fig. [2] and Fig.Fig. [6] with resemble the figures for and ki. 
The power-law correlation between s and k is reveal by Fig. HI where we fix 
m = p = 1 and tune 6 and n as we need. . 
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Fig. 7. Power-law probability distribution of weight w, the size of the network is 
5000 and the result is average of 10 individual experiments. Various values for p, 
m, n, 6 are chosen to display the interplay of different parameters. 

The evolution and distribution of weight can be calculated similarly as we deal 
with strength. Combine Eq. ([9]), Eq. (ITT!) ,Eq. (fT3|) . and define a = m^ + 2p5 + 
2n + 2p, we get 



dwij 2p6 + 2n w, 




dt a 
we can integrate the above equation and get 

m?+2pS + 2n 

Wij ~ tm2+2pa + 2n + 2p (^19) 

for large t. Therefore P{w) can be represented as P{w) = w~"'^ with 7^^, = 
2 + jy^2^25+2n - '^^^ simulations of Wij and P{w) are given representatives in 
Fig. i and Fig. [3 



5 Clustering Coefficients and Assortativeness 



Clustering coefficients depict connectivity among neighborhood of given ver- 
tices. The local clustering coefficient q for specific vertex i is defined to be 
Ci = ^ J2j<k sgn{wijWjkWik), which denotes the percentage of number of tri- 

angles in the neighborhood of i to the number of potential triangles there. Cj 
reveals how close vertices in the neighborhood of i are related. The average q 
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Fig. 8. (a) Clustering coefficients for different degree k. (b) Average nearest-neighbor 
degree for different degree k. The size of the network is 5000 and the result is average 
of 10 individual experiments. Various values for p, m, n, 6 are chosen to display the 
interplay of different parameters. 

of all vertices is denoted as C, and the average of Cj for vertices with degree k 
is denoted as C{k). Average degree of nearest neighbor knn,i for vertex i is 
also studied, as well as knn{k) for the average of knn, i of vertices with degree 
of k. knn reveals the assortativeness of a network. 

We perform numerical experiment the simulate the growth of network and 
analyze the clustering coefficient and assortativeness under the synthetical 
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weights' dynamic mechanism. In Fig. [8]^a) we show how clustering coefficients 
vary with tunable variables. The scale-free attributes C measures the overall 
density of triangles in the network. We can see that C is consistent with the 
network size while varying positively when those tunable variables change. 
The evolution of knn is show in Fig. [8]^b) which suggest the tunable assorta- 
tiveness of the network. 



6 Conclusion 

In summary, we propose a new model for weighted networks which synthesizes 
three sources of weights' dynamic: local weights' rearrangement raised by in- 
troduction of new vertex; self-increment of weights according to weights' pref- 
erential strategy; weights' creation and reinforcement proportional to strengths 
of both ends of nodes. Three sources independently contribute to the evolution 
and in the mean time also cooperatively interact. The homogenous behaviors 
that these weights' dynamics display suggest there may be some common un- 
derlying mechanisms that are not yet well understood. This model would be 
a good start for a synthetical and general understanding of weights' dynamic 
and hopefully our work would be helpful for the future study. 
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